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ABSTRACT 

Detrital  zircon  ages  from  quartzite  and 
metaconglomerate  in  the  Tusas  and  Picuris 
Mountains  in  northern  New  Mexico  reveal 
new  information  about  age  and  provenance 
trends  within  a  >1000  km2  Proterozoic  sedi¬ 
mentary  basin  and  provide  a  critical  test 
of  regional  correlations.  Samples  from  the 
Paleoproterozoic  Vadito  and  Hondo  groups 
are  dominated  by  a  single  detrital  zircon 
population  with  age  probability  peaks  that 
range  from  1765  to  1704  Ma.  Minor  Archean 
and  ca.  1850  Ma  age  probability  peaks  were 
also  recognized  in  some  samples.  Close  cor¬ 
respondence  between  detrital  zircon  ages 
and  the  age  of  surrounding  basement  rocks 
indicates  predominately  local  sources,  and 
we  interpret  systematic  shifts  in  peak  ages 
with  stratigraphic  position  to  represent 
changes  in  local  sources  through  time.  Simi¬ 
larities  of  age  spectra  support  previous  cor¬ 
relation  of  stratigraphic  units  between  dis¬ 
continuous  exposures  of  the  Hondo  Group. 
We  interpret  that  these  supracrustal  rocks 
were  deposited  in  a  single  basin  that  we  refer 
to  as  the  Pilar  basin. 

Two  samples  of  the  Marquenas  Forma¬ 
tion,  a  pebble  to  boulder  conglomerate  previ¬ 
ously  correlated  with  the  ca.  1700  Ma  Vadito 
Group,  are  dominated  by  Paleoproterozoic 
detrital  zircon  with  age  probability  peaks  at 
1707  and  1715  Ma  in  the  middle  and  upper 
units,  respectively.  Unlike  the  Vadito  and 
Hondo  group  samples,  the  Marquenas  For¬ 
mation  also  contains  abundant  ca.  1700- 
1600  Ma  zircon  derived  from  Mazatzal-aged 


sources  to  the  south  and  Mesoproterozoic  zir¬ 
con  with  age  probability  peaks  at  1479  and 
1457  Ma.  Weighted  averages  of  1477  ±  13  Ma 
and  1453  ±  10  Ma  for  the  youngest  detrital 
zircon  in  the  middle  and  upper  Marquenas 
Formation  provide  new  maximum  deposi- 
tional  age  constraints,  indicating  that  it  is 
not  part  of  the  Vadito  Group.  The  minimum 
age  is  not  well  constrained  but  is  interpreted 
to  be  ca.  1435  Ma  on  the  basis  of  the  timing 
of  regional  metamorphism  and  deformation 
previously  documented  in  the  Picuris  Moun¬ 
tains.  These  data  represent  the  first  evidence 
of  sedimentation  directly  associated  with  ca. 
1.4  Ga  regional  metamorphism,  plutonism, 
and  deformation  in  the  southwestern  United 
States  and  provide  an  important  new  con¬ 
straint  on  the  tectonic  evolution  of  southern 
Laurentia  during  this  time. 

INTRODUCTION 

Thick  metasedimentary  successions  that 
include  abundant  quartz  arenite  are  exposed 
throughout  the  southwestern  United  States 
and  in  many  Paleoproterozoic  orogens  around 
the  world  (e.g.,  Dott,  1983;  Rosen  et  al.,  1994; 
Cox  et  al.,  2002;  Medaris  et  al.,  2003;  Betts  and 
Giles,  2006;  Rainbird  and  Davis,  2007;  Jones 
et  al.,  2009).  Regional  and  global  correlation  of 
these  distinctive  sequences  provides  key  piercing 
points  for  proposed  Proterozoic  supercontinent 
reconstructions  (e.g.,  Dalziel,  1991;  Burrett  and 
Berry,  2000;  Sears  and  Price,  2000;  Karlstrom 
et  al.,  2001;  Goodge  et  al.,  2008)  and  constrains 
tectonic  models  for  the  growth  and  stabilization 
of  continental  lithosphere  (Giles  et  al.,  2002; 


Jones  et  al.,  2009).  In  the  southwestern  United 
States  (Fig.  1),  quartzite  and  related  rhyolite  suc¬ 
cessions  are  key  marker  units  for  distinguish¬ 
ing  the  age  and  tectonic  setting  of  events  that 
resulted  in  growth  of  southern  Laurentia  during 
the  Paleoproterozoic  (Karlstrom  and  Bowring, 
1988;  Hoffman,  1988;  Williams,  1991;  Jessup 
et  al.,  2006;  Whitmeyer  and  Karlstrom,  2007; 
Amato  et  al.,  2008;  Jones  et  al.,  2009). 

Perhaps  the  single  greatest  impediment  to 
understanding  the  tectonic  significance  of  these 
quartzite  successions  involves  uncertainties 
about  the  timing  of  deposition.  Absolute  ages 
are  only  available  for  a  small  number  of  locali¬ 
ties  (e.g.,  Bauer  and  Williams,  1989;  Cox  et  al., 
2002),  and  the  limited  geographical  extent  of 
individual  quartzite  exposures  complicates  the 
correlation  of  the  successions  at  the  regional 
and  global  scale.  Cross-cutting  relationships 
observed  throughout  the  southern  Rocky  Moun¬ 
tains  (Fig.  1)  indicate  that  deposition  occurred  in 
cycles  following  regional  orogenesis  at  ca.  1.70 
and  1.65  Ga  (Amato  et  al.,  2008;  Jones  et  al., 
2009).  However,  the  duration  of  sediment  depo¬ 
sition,  local  and  regional  basin  geometries,  and 
the  structural  and  metamorphic  evolution  of  the 
successions  are  not  well  understood. 

Here  we  present  new  detrital  zircon  geo¬ 
chronology  for  rocks  from  the  Paleoprotero¬ 
zoic  Vadito  and  Hondo  groups  in  northern  New 
Mexico  (Bauer  and  Williams,  1989;  Bauer, 
1993).  The  Vadito  and  Hondo  groups  consist 
of  >3  km  of  metasedimentary  and  metavolca- 
nic  rocks  that  are  exposed  in  multiple  mountain 
ranges  over  an  area  greater  than  1000  km2  (Figs. 
1  and  2).  These  rocks  are  broadly  correlated 
with  multiple  quartzite  successions  exposed  to 
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the  north  in  Colorado  (Jones  et  al.,  2009)  and  to 
the  west  in  Arizona  (Cox  et  al.,  2002),  and  their 
outcrop  trend  suggests  that  they  were  deposited 
in  an  elongated  depositional  basin  we  refer  to 
as  the  Pilar  basin  (Fig.  2).  We  interpret  the  Pilar 
basin  to  represent  one  of  multiple  northeast- 
striking,  broadly  contemporaneous  basins  that 
formed  in  southern  Colorado,  northern  New 
Mexico,  and  central  Arizona  ca.  1.7  Ga  follow¬ 
ing  the  culmination  of  the  Yavapai  orogeny  (Fig. 
1;  Jones  et  al.,  2009).  The  correlation  of  units 
within  the  Vadito  and  Hondo  groups  was  revised 
and  formalized  by  Bauer  and  Williams  (1989) 
and  is  widely  accepted  in  the  region.  However, 
provenance  studies  are  only  available  for  a  small 
number  of  localities,  and  correlations  within 
the  groups  are  complicated  by  uncertainty  in 
absolute  and  relative  age  constraints  and  poorly 
exposed  contact  relationships.  Our  new  results 
confirm  some  of  the  previous  correlations  and 
provide  new  insight  into  the  provenance  of  the 
thick  metasedimentary  succession.  However, 
our  findings  also  require  an  important  revision 
to  the  Proterozoic  lithostratigraphy  of  the  region 
and  may  also  have  direct  bearing  on  models  for 
regional  metamorphism  and  deformation  during 
the  Mesoproterozoic. 

GEOLOGICAL  SETTING 

Exposures  of  Precambrian  rocks  across  the 
southwestern  United  States  show  a  diverse 
assemblage  of  metavolcanic  and  metasedi¬ 
mentary  rocks  with  plutons  that  range  from 
mafic  to  felsic  composition.  These  rocks  were 
formed  and  accreted  to  the  southern  margin  of 
the  Archean  Wyoming  province  between  1800 
Ma  and  1600  Ma  (Condie,  1982;  Karlstrom  and 
Bowring,  1988;  Reed  et  al.,  1993)  as  part  of  a 
protracted  period  of  Laurentian  crustal  growth 
(Whitmeyer  and  Karlstrom,  2007).  These  expo¬ 
sures  have  been  divided  into  several  orogenic 
provinces  on  the  basis  of  rock  ages  and  isotopic 


Ligure  1.  Paleoproterozoic  quartzite  suc¬ 
cessions  of  the  southern  Rocky  Mountains, 
U.S.A.  Outcrop  localities  mentioned  in  the 
text  are  labeled  accordingly.  Faults  shown 
without  thrust  teeth  are  either  inferred 
thrust  faults  or  strike-slip  faults.  The  out¬ 
lines  of  possible  Paleoproterozoic  dep- 
ocenters  are  from  Jones  et  al.  (2009)  and 
highlight  ca.  1.70  Ga  metasedimentary  suc¬ 
cessions  that  might  be  broadly  correlative. 
The  Pilar  basin  discussed  in  the  text  is  at 
the  southern  end  of  the  southern  Colorado- 
northern  New  Mexico  depocenter  outline 
and  is  highlighted  in  Figure  2. 
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characteristics  (Fig.  1).  The  Yavapai  province 
is  interpreted  to  represent  a  complex  collage  of 
predominantly  juvenile  arc  terranes  character¬ 
ized  by  rocks  with  Nd  model  ages  between  2000 
and  1800  Ma  (Bennett  and  DePaolo,  1987). 
Rocks  of  the  Yavapai  province  were  accreted 
to  the  Laurentian  margin  between  1780  and 
1700  Ma  along  a  belt  stretching  from  Colorado 
to  Arizona  and  New  Mexico.  The  final  colli- 


sional  phase  of  this  long-lived,  progressive  oro- 
genic  event  occurred  ca.  1710-1700  Ma  (Karl- 
strom  and  Bowring,  1988),  and  was  followed 
by  voluminous  postorogenic  granitoid  magma- 
tism  from  ca.  1690  to  1650  Ma  (Anderson  and 
Cullers,  1999). 

The  Mazatzal  province  lies  south  of  the 
Yavapai  province  and  extends  across  central 
and  southern  New  Mexico  and  Arizona  (Fig.  1). 


Mazatzal  province  rocks  are  characterized  by  Nd 
model  ages  between  1800  and  1700  Ma  (Bennett 
and  DePaolo,  1987)  and  were  accreted  to  south¬ 
ern  Laurentia  during  the  1660-1600  Ma  Mazat¬ 
zal  orogeny  (Silver,  1965;  Karlstrom  and  Bow¬ 
ling,  1988;  Luther,  2006;  Amato  et  al.,  2008). 
Mazatzal-aged  deformation  affected  a  large 
foreland  region  of  the  southern  Yavapai  province 
(Transition  Zone,  Fig.  1;  Karlstrom  and  Hum¬ 
phreys,  1998),  and  the  Mazatzal  Deformation 
Front  represents  the  approximate  northern  extent 
of  these  effects  (Shaw  and  Karlstrom,  1999). 
Various  workers  have  challenged  the  juvenile 
arc  accretion  model  for  the  Yavapai  and  Mazat¬ 
zal  orogenies  on  the  basis  of  zircon  ages,  litho¬ 
logical  associations,  and  limited  Hf  isotopic  data 
(Bickford  and  Hill,  2007a;  Bickford  et  al.,  2008), 
but  alternative  models  are  still  being  evaluated 
and  debated  (Duebendorfer,  2007;  Karlstrom  et 
al.,  2007;  Bickford  and  Hill,  2007b). 

After  a  ca.  100  m.y.  tectonic  lull  from  ca. 
1600  to  1500  Ma,  renewed  southward  growth  of 
Laurentia  is  inferred  to  have  occurred  during  the 
Mesoproterozoic.  This  interpretation  is  based 
upon  Nd  model  ages  of  1500-1300  Ma  from 
the  granite-rhyolite  province  that  extends  from 
northern  Mexico  to  Labrador,  Canada  (Bennett 
and  DePaolo,  1987;  Patchett  and  Ruiz,  1989; 
Karlstrom  et  al.,  2001).  An  episode  of  wide¬ 
spread  granitic  magmatism  and  local  emplace¬ 
ment  of  mafic  dikes  occurred  throughout  the 
southwestern  United  States  between  1470  and 
1360  Ma  (Williams,  1991;  Reed  et  al.,  1993; 
Williams  et  al.,  1999),  and  rocks  of  this  age  cur¬ 
rently  account  for  nearly  20%  of  all  Precambrian 
exposures  across  the  region.  High  crustal  tem¬ 
peratures  of  300-750  °C  during  the  same  time 
caused  regional  greenschist-  to  near  granulite- 
facies  metamorphism  (Grambling  et  al.,  1989; 
Williams  and  Karlstrom,  1996;  Williams  et  al., 
1999;  Daniel  and  Pyle,  2006)  and  reset  argon 
ages  throughout  the  region  (Grambling  and 
Dallmeyer,  1993;  Karlstrom  et  al.,  1997;  Shaw 
et  al.,  1999,  2005;  Sanders  et  al.,  2006).  Circa 
1.4  Ga  granites,  previously  described  as  being 
A-type  because  of  their  alkalinity,  anhydrous 
character,  and  presumed  anorogenic  tectonic 
setting  (Loiselle  and  Wones,  1979;  Anderson, 
1983),  are  ferroan  in  nature  (Frost  et  al.,  2001), 
a  geochemical  characteristic  that  is  indicative 
of  mantle  influence  (Frost  and  Frost,  1997)  and 
is  generally  associated  with  extensional  tec¬ 
tonic  environments  such  as  continental  rifting 
(Emslie,  1978;  Whalen  et  al.,  1987;  Eby,  1990). 

However,  regional  evidence  exists  for  short¬ 
ening  or  strike-slip  deformation  within  the  ther¬ 
mal  aureoles  of  plutons  and  contemporaneous 
reactivation  of  northeast-striking  crustal  shear 
zones  in  the  Rocky  Mountains  and  southwestern 
United  States  (Graubard  and  Mattinson,  1990; 


Figure  2.  Reconstruction  of  Proterozoic  orogenic  belt  across  north-central  New  Mexico 
(after  Karlstrom  and  Daniel,  1993;  Cather  et  al.,  2006).  Exposures  of  Hondo  and  Vadito 
Group  rocks  correspond  to  regional  aeromagnetic  lows  that  we  interpret  to  represent  the 
extent  of  the  Pilar  basin  development.  Outcrop  localities  discussed  in  text  are  labeled  on  the 
inset  and  map  and  reconstruction:  C — Cimarron  Mountains;  N — Nacimiento  Mountains; 
Pic — Picuris  Mountains;  R — Rincon  Mountains;  RM — Rio  Mora;  SF — Santa  Fe  Range; 
Ta — Taos  Range;  Tr — Truchas  Peaks;  TU — Tusas  Mountains. 
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Shaw  et  al.,  2001;  McCoy  et  al.  2005;  Jessup  et 
al.,  2006;  Jones  et  al.,  2010).  Kinematic  evidence 
for  regional,  northwest-directed  crustal  shorten¬ 
ing  synchronous  with  pluton  emplacement  and 
regional  metamorphism  suggests  that  ca.  1.4  Ga 
magmatism  coincided  with  regional  orogenesis 
associated  with  a  convergent  plate  boundary 
along  the  distal  southern  margin  of  Laurentia 
(e.g.,  Nyman  et  al.,  1994).  These  events  were 
all  part  of  a  prolonged,  ca.  800  m.y.  episode  of 
crustal  growth  along  the  long-lived  “southern” 
margin  of  Laurentia  that  culminated  in  the  Gren¬ 
ville  orogeny  and  assembly  of  the  supercontinent 
Rodinia  at  ca.  1.1  Ga  (Karlstrom  et  al.,  2001). 

Proterozoic  Quartzite  Successions 

During  the  southward  growth  of  Laurentia, 
successions  of  quartz  sandstone  up  to  2  km 
thick  were  deposited  within  the  orogenic  belts. 
Quartzite  successions  occur  in  both  the  Yavapai 
and  Mazatzal  provinces,  and  include  exten¬ 
sively  exposed  units  like  the  Ortega  Formation 
(New  Mexico)  and  Uncompahgre  Formation 
(Colorado)  as  well  as  numerous  smaller,  more 
localized  exposures  with  similar  lithologies  and 
outcrop  characteristics  (Fig.  1).  Existing  age 
constraints  suggest  that  sedimentation  occurred 
ca.  1700-1680  Ma  following  the  Yavapai  orog¬ 
eny  (Jones  et  al.,  2009)  and  ca.  1650-1600  Ma 
following  the  Mazatzal  orogeny  (Luther  et  al., 
2005a,  2006;  Luther,  2006;  Amato  et  al.,  2008). 
Quartzite  is  commonly  nearly  pure  (>95% 
quartz)  with  minor  muscovite,  Al-silicates, 
hematite,  zircon,  and  monazite.  Primary  sedi¬ 
mentary  structures  are  locally  well  preserved, 
including  common  cross  stratification.  Deposi- 
tional  facies  are  similar  from  bottom  to  top  and 
region  to  region  and  indicate  shallow  marine 
(<10  m  water  depth)  or  fluvial  environments 
(Trevena,  1979;  Soegaard  and  Eriksson,  1985, 
1989;  Harris  and  Eriksson,  1990). 

In  many  outcrop  localities,  quartzite  directly 
overlies  thick  successions  of  voluminous,  high- 
silica  rhyolite  (Fig.  1).  The  contact  between 
rhyolite  and  quartzite  is  generally  interlayered 
to  gradational  and  is  commonly  marked  by  a 
distinctive,  Mn-rich  contact  interval  (Bauer 
and  Williams,  1989).  Basement  rock  assem¬ 
blages  underlying  quartzite-rhyolite  successions 
regionally  are  typically  characterized  by  meta¬ 
morphosed  mafic  volcanic  rocks  and  volcano- 
genic  marine  metasedimentary  rocks  (Bauer  and 
Williams,  1989;  Jessup  et  al.,  2005).  Examples 
of  these  older  (1800-1720  Ma;  Condie,  1982) 
assemblages  include  the  Pecos,  Gold  Hill,  and 
Moppin  complexes  in  northern  New  Mexico 
(Bauer  and  Williams,  1989)  and  the  Dubois  and 
Cochetopa  successions  in  southern  Colorado 
(Bickford  and  Boardman,  1984).  Basement 


rocks  commonly  contain  evidence  for  multiple 
episodes  of  deformation  and/or  metamorphism 
that  are  not  recognized  in  the  overlying  quartz¬ 
ite  (Gibson  and  Harris,  1992).  Thus,  the  contact 
between  basement  assemblages  and  quartzite 
successions  has  been  variably  interpreted  as  an 
unconformable  depositional  contact,  a  sheared 
unconformity,  or  a  fault  contact. 

The  quartzite-rhyolite  successions  and  under¬ 
lying  basement  assemblages  were  deformed  by 
folding  and  thrust  imbrication.  Across  Colorado, 
quartzite  exposures  are  found  as  tight,  upright 
synclinal  “keels”  interpreted  to  be  the  roots 
of  larger,  now  eroded  folds  (Wells  et  al.,  1964; 
Reuss,  1974;  Jones  et  al.,  2009).  In  the  Needle 
Mountains  of  southwestern  Colorado  and  Tusas 
Mountains  of  northern  New  Mexico,  1-2-km- 
thick  sections  of  quartzite  define  tight  to  open, 
large-wavelength  (kilometer-scale)  folds  con¬ 
sistent  with  some  fold-and-thrust-belt  geom¬ 
etries  (Harris,  1990;  Williams,  1991;  Williams 
et  al.,  1999).  Deformation  and  metamorphism 
of  quartzite  successions  in  the  southern  Yavapai 
province  are  interpreted  to  have  occurred  during 
accretion  of  the  Mazatzal  province  to  the  south 
(Jones  et  al.,  2009).  The  metasedimentary  suc¬ 
cessions  were  metamorphosed  at  greenschist- 
to  amphibolite-facies  conditions  prior  to  ca. 
1450  Ma,  when  they  were  widely  intruded  by 
coarse-grained  granitic  plutons  (Barker,  1969; 
Reuss,  1974;  Jessup  et  al.,  2006;  Jones  and  Con¬ 
nelly,  2006).  Deformation  and  metamorphism  of 
Mazatzal  province  successions  possibly  occurred 
during  the  waning  stages  of  the  Mazatzal  orog¬ 
eny  (Luther  et  al.,  2005a;  Amato  et  al.,  2008)  but 
may  have  also  occurred  during  the  Mesoprotero¬ 
zoic  (Baer  et  al.,  2003;  Amato  et  al.,  201 1). 

PROTEROZOIC  GEOLOGY  OF 
NORTHERN  NEW  MEXICO 

Lithostratigraphy 

Exposures  of  Proterozoic  rocks  in  north¬ 
ern  New  Mexico  are  dominated  by  granite- 
greenstone  associations,  supracrustal  succes¬ 
sions,  and  granitoids  (Karlstrom  et  al.,  2004). 
The  oldest  units  include  the  ca.  1770-1750  Ma 
Moppin  and  Gold  Hill  complexes  exposed  in  the 
Tusas  Mountains  and  Taos  Range,  respectively 
(Figs.  1  and  2;  Bauer  and  Williams,  1989,  and 
references  therein),  and  the  ca.  1720  Ma  Pecos 
complex,  exposed  in  the  eastern  Santa  Fe  Range 
(Fig.  2).  These  complexes  are  dominated  by 
mafic  metavolcanic  rocks  but  also  include  fel- 
sic  metavolcanic  rocks,  metasedimentary  rocks, 
and  iron  formation  (Bauer  and  Williams,  1989; 
Robertson  and  Condie,  1989;  Karlstrom  et  al., 
2004).  They  are  represented  in  black  in  Figure 
2,  and  they  are  widely  intruded  by  calc-alkaline 


tonalite  to  granodiorite  plutons  (purple  in  Fig. 
2).  The  Moppin  and  Gold  Hill  complexes  are 
interpreted  to  have  been  formed  in  arcs  in  a 
marine  setting  (Condie,  1986).  Robertson  and 
Condie  (1989)  interpreted  that  the  Pecos  com¬ 
plex  was  formed  in  a  backarc  setting,  possibly 
in  a  continental  margin  setting. 

The  Vadito  and  Hondo  groups  are  successions 
of  metavolcanic  and  metasedimentary  rocks 
that  are  in  unconformable  or  tectonic  contact 
with  the  older  mafic  complexes.  These  succes¬ 
sions  are  exposed  in  mountain  ranges  through¬ 
out  northern  New  Mexico  and  are  represented 
in  yellow  and  green  in  Figures  2  and  3.  Units 
within  the  Vadito  Group  are  not  well  defined 
across  northern  New  Mexico  (Fig.  4;  Bauer  and 
Williams,  1989).  Stratigraphic  relationships  are 
obscured  by  discontinuous  exposure,  a  seem¬ 
ingly  greater  degree  of  stratigraphic  variation 
along  strike  (Bauer  and  Williams,  1989),  and 
multiple  phases  of  deformation  that  disrupted 
and  possibly  repeated  parts  of  the  original  stra¬ 
tigraphy  (Fig.  3).  In  the  Picuris  Mountains,  the 
lowermost  section  of  the  Vadito  Group  largely 
consists  of  amphibolite  with  minor  quartz- 
biotite-muscovite  schist,  cross-bedded  impure 
quartzite,  and  metaconglomerate.  To  the  north, 
various  schist  and  minor  quartzite  are  dominant. 
Rare  cross  stratification  within  quartzite  lenses 
indicates  stratigraphic  younging  to  the  north. 
McCarty  (1983)  interpreted  these  rocks  as 
interlayered  graywacke,  quartzite,  volcaniclas- 
tic  sediment,  and  basaltic  flows,  and  inferred  a 
deep-water  basin  setting  for  deposition. 

The  500-m-thick  Marquenas  Formation 
is  exposed  to  the  north  of  schists  of  the  lower 
Vadito  Group  in  the  Picuris  Mountains  (Figs.  3 
and  4B).  It  consists  of  a  lower,  poorly  sorted, 
highly  strained,  polymictic  boulder  conglom¬ 
erate,  a  middle  quartzite,  and  an  upper,  highly 
strained,  pebble-to-cobble  conglomerate.  Cross 
stratification  within  the  Marquenas  Formation 
indicates  that  the  layering  is  overturned  and  that 
the  unit  youngs  to  the  north,  supporting  strati¬ 
graphic  interpretations  that  place  the  Marque¬ 
nas  Formation  within  the  Vadito  Group  (Fig. 
4B;  Bauer  and  Williams,  1989).  Previous  work 
(Miller  et  al.,  1963;  Grambling,  1979;  Gram- 
bling  and  Codding,  1982)  in  the  Truchas  Peaks 
and  Rio  Mora  areas  (Fig.  2)  east  of  the  Picuris 
Mountains  documented  a  relatively  thick 
sequence  of  conglomerate  that  is  correlated  with 
the  Marquenas  Formation. 

The  uppermost  part  of  the  Vadito  Group  is 
represented  by  the  Glenwoody  Formation  in  the 
Picuris  Mountains  (Fig.  4B).  This  unit  comprises 
-300  m  of  feldspathic  schist  that  is  exposed 
beneath  the  Ortega  Formation  in  the  northern 
Picuris  Mountains  (Bauer  and  Williams,  1989). 
Despite  poor  exposure  of  the  upper  and  lower 
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contact,  Bauer  and  Williams  (1989)  interpreted 
the  Glenwoody  Formation  to  represent  the  upper 
Vadito  Group  (Fig.  4B)  locally  on  the  basis  of 
similarities  with  other  upper  Vadito  Group  rocks 
in  other  ranges  and  the  presence  of  a  Mn-rich 
horizon  that  is  typically  found  high  in  the  Vadito 
Group  stratigraphy.  Bauer  et  al.  (2005)  mapped 
transitional  rocks  between  the  Vadito  and  Hondo 
groups  in  the  eastern  Picuris  Mountains.  These 
transitional  units  are  described  as  feldspathic 
schists  of  the  Vadito  Group  that  grade  later¬ 
ally  into  various  quartzite,  conglomerate,  and 
schistose  quartzite  and  conglomerate  units  that 
grade  upward  into  the  massive  orthoquartzite  of 
the  Ortega  Formation  (Bauer  et  al.,  2005).  The 
stratigraphy  here  is  disrupted  by  the  Plomo  fault, 
but  Bauer  et  al.  (2005)  suggested  that  these  rocks 
are  equivalent  to  similar  transitional  rocks  in  the 
Tusas  Mountains  (Fig.  4). 

In  the  Tusas  Mountains  to  the  west  (Figs.  1  and 
2),  Vadito  Group  exposures  are  dominated  by 
immature  quartzofeldspathic  metasedimentary 
rocks  and  felsic  metavolcanic  rocks  (Bauer  and 
Williams,  1989).  These  rocks  range  from  highly 
schistose  to  massive  and  are  interlayered  with 
amphibolite,  cross-bedded  muscovite  quartzite, 
feldspathic  quartzite,  biotite-muscovite  schist, 
conglomerate,  and  pelitic  schist  (Bauer  and 
Williams,  1989,  and  references  therein).  The 
contacts  between  these  different  rock  types  are 


Figure  3.  Simplified  geologic  map  of  the 
Picuris  Mountains  and  cross  section  (A-A') 
across  the  western  part  of  the  range  showing 
the  relationship  of  the  Marquenas  Forma¬ 
tion  to  the  Vadito  and  Hondo  groups  (modi¬ 
fied  from  Bauer,  1993;  Bauer  and  Helper, 
1994).  Approximate  sample  locations  of 

(1)  J07-PIC1,  (2)  J07-PIC2,  (3)  J07-PIC3, 
(4)  J07-PIC4,  and  (5)  J07-PIC5  are  shown 
in  the  cross  section  in  the  appropriate  struc¬ 
tural  position.  Summary  stereonet  plots  for 
the  area  surrounding  the  sampled  Marque¬ 
nas  Formation  exposures  show  (A)  average 
orientation  of  SO  and  S2  in  the  Marquenas 
Formation  after  McCarty  (1983);  (B)  aver¬ 
age  orientation  of  SO  and  S2  in  the  Vadito 
and  Hondo  groups  after  Holcombe  and  Cal- 
lendar  (1982);  (C)  average  orientation  of 
SO  and  S2  in  the  Vadito  and  Hondo  groups 
after  McCarty  (1983);  and  (D)  average 
orientation  of  S0/S1  and  S2  in  the  Hondo 
Group  together  with  the  average  LI  stretch¬ 
ing  lineation  on  S0/S1  in  the  Hondo  Group 
after  Bauer  (1988).  The  contours  show  the 
X-axis  (stretching  lineation)  of  clasts  in  the 
Marquenas  Formation  after  Nielsen  and 
Scott  (1979). 


poorly  exposed  and  are  obscured  by  polyphase 
deformation  that  involved  isoclinal  folding 
and  development  of  penetrative  metamorphic 
foliations.  One  notable  unit  within  the  Vadito 
Group  stratigraphy  in  the  Tusas  Mountains  is 
the  ~60-m-thick  Big  Rock  Formation  (Fig.  4A). 
This  unit  is  divided  into  a  lower  member  domi¬ 
nated  by  conglomerate  with  minor  interlayered 
micaceous  quartzite  and  an  upper  member  dom¬ 
inated  by  micaceous  and  feldspathic  quartzite 
(Bauer  and  Williams,  1989).  Metaconglomerate 
layers  contain  elongate  cobbles  up  to  20  cm  in 
diameter,  and  the  cobble  compositions  include 
vein  quartz,  quartzite,  schist,  and  metarhyolite. 
The  clast  size  is  smaller  on  average  than  the 
Marquenas  Formation,  and  the  Marquenas  For¬ 
mation  metaconglomerate  is  thicker  and  is  inter¬ 


layered  with  a  larger  volume  of  impure  quartz¬ 
ite.  Otherwise,  the  two  units  are  lithologically 
similar  and  display  younging  criteria  such  as 
cross  beds  and  graded  bedding  suggesting  that 
they  are  stratigraphically  beneath  the  Hondo 
Group.  Thus,  they  were  tentatively  correlated 
within  the  upper  Vadito  Group  by  Bauer  and 
Williams  (1989)  and  Bauer  (1993)  (Fig.  4). 

The  upper  Vadito  Group  in  the  Tusas  Moun¬ 
tains  contains  the  Burned  Mountain  Formation 
(Fig.  4A).  It  is  dominated  by  metarhyolite  that 
contains  distinctive  quartz  and  feldspar  eyes  in  a 
fine-grained,  laminated  matrix  (Bauer  and  Wil¬ 
liams,  1989).  The  formation  ranges  in  thickness 
from  4.5  to  30  m  and  yielded  a  U-Pb  zircon  age 
ofca.  1700  Ma  (Bauer  and  Wilhams,  1989).  This 
age  is  interpreted  to  represent  crystallization 
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Figure  4.  Interpretive  structural  and  stratigraphic  reconstructions  of  Proterozoic  lithostrati- 
graphic  sequences  in  the  Tusas  and  Picuris  Mountain,  New  Mexico.  (A)  Lithostratigraphic 
relationships  of  the  Vadito  and  Hondo  groups,  Tusas  Mountains,  northern  New  Mexico 
(Bauer  and  Williams,  1989;  Williams,  1991).  (B)  Lithostratigraphy  of  the  Vadito  and  Hondo 
groups,  Picuris  Mountains,  northern  New  Mexico  (modified  from  Bauer,  1993).  (C)  Rein¬ 
terpreted  lithostratigraphy  of  (B)  showing  the  Mesoproterozoic  Marquenas  Formation  in 
unconformable  depositional  contact  or  fault  contact  with  the  Paleoproterozoic  Vadito  and 
Hondo  groups. 
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of  the  rhyolitic  protolith,  thus  providing  one 
of  the  only  constraints  on  the  absolute  age  of 
the  Vadito  and  Hondo.  The  uppermost  Vadito 
Group  in  the  Tusas  Mountains  is  characterized 
by  immature  micaceous  quartzite  (Transitional 
rocks  in  Fig.  4A)  that  grades  upward  into  more 
thickly  bedded  orthoquartzite  of  the  lower 
Ortega  Formation. 

The  Hondo  Group  (Fig.  4)  overlies  the  Vadito 
Group  and  consists  of  several  metasedimentary 
formations.  The  basal  Ortega  Formation,  an 
~l-km-thick,  cross-bedded  quartzite,  is  overlain 
by  interlayered  pelitic  schist  and  cross-bedded 
quartzite  of  the  Rinconada  Formation.  Upsec- 
tion,  the  Pilar  Formation  consists  of  a  black 
graphitic  phyllite  followed  by  graphitic,  gamet- 
biotite  schist,  phyllite  and  micaceous  quartzite 
of  the  Piedra  Lumbre  Formation.  Soegaard  and 
Eriksson  (1986)  interpreted  the  Hondo  Group  to 
represent  fluvial  to  shallow-marine  deposition 
on  a  southward-deepening  siliciclastic  shelf. 

Proposed  reconstructions  of  right-slip  faulting 
across  the  Tusas,  Picuris,  Truchas,  and  Rio  Mora 
ranges  (Karlstrom  and  Daniel,  1993;  Daniel  et 
al.,  1995;  Cather  et  al.,  2006)  bring  Protero¬ 
zoic  rock  units  and  deformational  features  into 
alignment  and  help  define  an  east-west  trend  of 
Vadito  and  Hondo  group  exposures  that  we  refer 
to  as  the  Pilar  basin  (Fig.  2).  The  Pilar  basin  sits 
just  north  of  a  pronounced  geophysical  linea¬ 
ment  interpreted  to  represent  a  major  Protero¬ 
zoic  accretionary  boundary  in  the  southwestern 
United  States  (Fig.  1 ;  Karlstrom  and  Humphreys, 
1998).  A  zone  of  relatively  weak  magnetism 
located  north  of  the  Chama  lineament  corre¬ 
sponds  well  with  the  surface  exposures  of  the 
Vadito  and  Hondo  groups  (Cather  et  al.,  2006). 
We  interpret  this  magnetic  low  to  represent  the 
continuation  of  these  groups  in  the  subsurface 
and,  therefore,  place  the  approximate  northern 
boundary  of  the  basin  along  its  northern  edge. 
Quartzite  exposures  north  of  this  boundary  in  the 
Cimarron  and  Taos  ranges  may  be  correlative 
with  the  Hondo  Group  (Fig.  2).  Furthermore,  the 
northwestern  trend  of  the  Pilar  basin  in  Figure 
2  suggests  that  the  Hondo  Group  might  also  be 
correlated  with  the  Uncompahgre  Formation  in 
the  Needle  Mountains  of  southwestern  Colorado 
(Fig.  1).  However,  these  correlations  have  not 
yet  been  tested,  and  the  quartzite-rhyolite  asso¬ 
ciation  of  the  Vadito  and  Hondo  groups  has  not 
been  recognized  in  these  areas. 

Deformation  and  Metamorphism  of 
the  Hondo  Group,  Vadito  Group,  and 
Marquenas  Formation 

Three  generations  of  ductile  deformation  are 
recognized  in  exposures  of  the  Hondo  Group, 
Vadito  Group,  and  Marquenas  Formation  in 


northern  New  Mexico  (Williams,  1991;  Bauer, 
1993;  Williams  et  al.,  1999).  D1  includes  local 
bedding-parallel  foliation  (SI)  and  fold  nappes 
(Williams,  1991;  Read  et  al.,  1999).  D2  pro¬ 
duced  the  dominant,  penetrative  foliation  (S2), 
north-verging,  km-scale  inclined  folds  (F2)  and 
ductile  faults  in  the  Picuris  and  Tusas  Moun¬ 
tains  (Nielsen  and  Scott,  1979;  Holcombe  and 
Callender,  1982;  McCarty,  1983;  Bauer,  1988, 
1993;  Williams,  1991).  D3  structures  include 
weakly  developed  upright,  open  folds  and  an 
associated  crenulation  cleavage  (S3),  but  the 
distribution  and  intensity  of  D3  deformation  is 
heterogeneous  throughout  the  region  (Williams 
et  al.,  1999).  Deformation  was  accompanied 
by  multiple  metamorphic  events  (M1-M3), 
culminating  with  regional  amphibolite-facies 
metamorphism  to  550-700  °C  ca.  1.47- 
1.42  Ga  (Williams  et  al.,  1999).  Map  relation¬ 
ships  around  dated  plutons  and  textural  studies 
of  dated  minerals  indicate  that  D3  deformation 
occurred  ca.  1.42  Ga  (Read  et  al.,  1999;  Wil¬ 
liams  et  al.,  1999;  Daniel  and  Pyle,  2006).  The 
ages  of  D1  and  D2  deformation  are  more  uncer¬ 
tain  and  could  be  ca.  1.65  or  1.42  Ga  (Williams 
et  al.,  1999;  Daniel  and  Pyle,  2006).  Williams 
(1991)  and  Bauer  (1993)  proposed  that  D2  rep¬ 
resents  a  progressive  north-vergent  deforma¬ 
tional  event  that  could  possibly  be  a  continua¬ 
tion  of  a  more  cryptic  D1  regional  deformation. 

In  the  Picuris  Mountains,  relict  bedding  (SO) 
in  the  Marquenas  Formation  strikes  east  and 
dips  steeply  south  (Fig.  3A;  McCarty,  1983), 
roughly  parallel  to  SO  in  Vadito  and  Hondo 
Group  rocks  to  the  north  and  south  (Figs.  3B  and 
3C;  Holcombe  and  Callender,  1982;  McCarty, 
1983).  Consistent  stratigraphic-facing  direction 
throughout  the  Marquenas  Formation  indicates 
that  the  unit  is  overturned  to  the  north  but  is 
not  internally  folded.  The  Marquenas  Forma¬ 
tion  does  not  contain  the  same  bedding-parallel 
schistosity  (SI)  that  is  observed  in  schistose 
units  of  the  Vadito  Group,  so  it  is  possible  that 
the  Vadito  and  Hondo  Groups  were  deformed 
and  metamorphosed  at  least  once  prior  to  F2 
folding.  But  the  Marquenas  Formation  does  con¬ 
tain  an  east-northeast-striking,  steeply  south¬ 
dipping  penetrative  tectonic  foliation  that  is 
similar  in  style  and  orientation  to  the  dominant, 
penetrative  S2  foliation  in  Vadito  and  Hondo 
Group  rocks  (Figs.  3A-3C;  Bauer,  1988, 1993). 
The  Plomo  fault  is  a  steeply  south-dipping  zone 
of  high  strain  that  separates  the  Marquenas 
Formation  from  folded  Hondo  Group  rocks  to 
the  north  (Bauer,  1993).  Quartzite  clasts  in  the 
Marquenas  Formation  are  strongly  flattened  in 
a  mylonitic  foliation  that  parallels  the  fault  con¬ 
tact,  and  elongated  clasts  define  a  prominent 
south-plunging  stretching  lineation  (Fig.  3D) 
with  reverse-sense  (top-up-to-the-north)  kine¬ 
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matics  (Bauer,  1993).  These  relationships  sug¬ 
gest  that,  at  a  minimum,  the  Marquenas  Forma¬ 
tion  was  deformed  together  with  the  Vadito  and 
Hondo  groups  during  megascopic  F2  folding 
and  ductile  shearing  along  the  Plomo  fault  in  the 
Picuris  Mountains.  The  shear  zone  cuts  across 
the  S2  foliation  and  north-vergent  F2  folds  (Fig. 
3;  Bauer,  1993;  Bauer  and  Helper,  1994),  indi¬ 
cating  that  some  ductile  movement  on  the  fault 
postdated  D2.  The  total  stratigraphic  separation 
across  the  Plomo  fault  is  not  well  established  but 
is  estimated  to  be  as  much  as  several  kilometers 
(Bauer,  1993). 

Metamorphic  mineral  assemblages  in  expo¬ 
sures  throughout  the  Picuris  Mountains  indicate 
amphibolite-facies  conditions  during  the  domi¬ 
nant  D2  deformation,  and  the  Marquenas  Forma¬ 
tion  locally  contains  garnet,  biotite,  and  chlorite 
in  thin,  schistose  partings.  Peak  metamorphic 
conditions  throughout  the  range  were  near 
the  aluminosilicate  triple  point  (500-550  °C, 
and  4.0  kb;  Holdaway,  1978;  McCarty,  1983; 
Grambling  and  Williams,  1985;  Read  et 
al.,  1999;  Williams  et  al.,  1999;  Daniel  and 
Pyle,  2006)  at  or  near  the  time  of  major  fold¬ 
ing  (Nielsen  and  Scott,  1979;  McCarty,  1983; 
Bauer,  1988,  1993).  Kyanite  and  sillimanite  are 
typically  deformed  and  are  thus  interpreted  to 
be  pre-  to  syntectonic  with  respect  to  F2  folding 
and  shearing  (Bauer,  1988,  1993).  Andalusite 
crystallized  after  most  of  the  folding,  indicating 
that  metamorphism  outlasted  deformation  in 
the  Picuris  Mountains  and  surrounding  ranges 
(Grambling  and  Williams,  1985;  Grambling  et 
al.,  1989;  Bauer,  1993;  Williams,  1991;  Bauer, 
1993).  Subhorizontal  metamorphic  isograds 
appear  to  overprint  map-scale  ductile  structures 
(Grambling,  1979;  Grambling  and  Williams, 
1985;  Bauer,  1993).  In  exposures  of  Vadito 
Group  schists  -50-100  m  south  of  the  Plomo 
fault  and  the  Marquenas  Formation,  Williams  et 
al.  (1999)  reported  metamorphic  conditions  of 
600-650  °C  and  3  kb  for  andalusite-  and  cordi- 
erite-bearing  rocks.  Monazite  inclusions  within 
andalusite  porphyroblasts  were  reported  to  be 
ca.  1450  Ma  (Williams  et  al.,  1999).  In  Hondo 
Group  rocks  exposed  in  the  northern  Picuris 
Mountains,  Daniel  and  Pyle  (2006)  documented 
peak  metamorphic  conditions  of  525-540  °C 
and  4.0-4.2  kb  that  occurred  between  1450  and 
1435  Ma.  The  disparity  in  pressures  and  tem¬ 
peratures  from  south  to  north  may  reflect  late 
movement  across  the  Plomo  fault.  Alternatively, 
the  pressure  gradient  may  reflect  a  difference 
in  erosion  level,  with  greater  uplift  and  erosion 
along  the  northern  margin  of  the  Picuris  Moun¬ 
tains  relative  to  the  south,  possibly  related  to 
uplift  along  the  Embudo  fault  zone.  Temperature 
differences  between  the  northern  and  southern 
Picuris  Mountains  may  partly  reflect  differences 
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in  the  petrogenetic  grids  used  to  estimate  tem¬ 
peratures.  Williams  et  al.  (1999)  used  Pattison’s 
(1992)  triple  point,  -550  °C  and  4.5  kbar,  versus 
the  triple  point  of  Holdaway  (1971),  -500  °C 
and  4  kbar,  used  by  Daniel  and  Pyle  (2006). 

DETRITAL  ZIRCON 
GEOCHRONOLOGY 

Methods 

We  present  new  U-Pb  detrital  zircon  age  data 
for  quartzite  and  conglomerate  collected  from 
exposures  in  the  Tusas  Mountains  and  Picuris 
Mountains  in  northern  New  Mexico.  Sample 
processing  followed  the  methodology  of  Jones 
and  Connelly  (2006).  Following  an  initial  mag¬ 
netic  separation  to  remove  paramagnetic  min¬ 
erals  and  zircon  with  abundant  inclusions,  the 
remaining  zircon  grains  were  handpicked  to 
include  populations  representing  each  of  the 
various  morphologies,  sizes,  and  colors  that 
were  recognized.  Approximately  150  grains 
from  each  sample  were  mounted  in  1-inch  epoxy 
discs  and  ground  and  lightly  polished  to  reveal 
the  grain  surfaces.  We  analyzed  zircon  at  the 
Geological  Survey  of  Denmark  and  Greenland 
(GEUS)  using  single-collector,  laser  ablation- 
magnetic  sector  field-inductively  coupled 
plasma-mass  spectrometry  (LA-SF-ICP-MS) 
techniques  described  by  Frei  and  Gerdes  (2009). 

Age  data  are  summarized  in  Table  1,  and  full 
analytical  data  and  sample  location  coordinates 
are  available  in  Supplemental  Table  l1.  Figures 
5  and  6  show  age-distribution  curves  and  U-Pb 
concordia  diagrams  for  all  detrital  zircon  from 
all  samples.  Figure  7  shows  age-distribution 
curves  of  <10%  discordant  Proterozoic  detrital 
zircon  for  all  samples  together  with  published 
detrital  zircon  data  from  Ortega  Formation 
quartzite  exposed  in  the  Tusas  Mountains,  New 
Mexico  (Fig.  7E;  Jones  et  al.,  2009)  and  zircon 
crystallization  ages  of  basement  rocks  in  central 
and  northern  New  Mexico  (Fig.  71;  Karlstrom 
et  al.,  2004).  Peak  ages  in  Figures  5-7  represent 
the  weighted  average  of  clusters  of  three  grains 
or  more  with  overlapping  ages  (Table  1). 

We  used  Kolmogorov-Smimoff  (K-S)  sta¬ 
tistics  (Press  et  al.,  1986)  to  help  assess  the 
similarity  of  the  detrital  zircon  age  populations 
among  the  samples  (Berry  et  al.,  2001 ;  Degraaff- 
Surpless  et  al.,  2003;  Dickinson  and  Gehrels, 
2009).  We  also  included  published  detrital 


'Supplemental  Table  1.  Excel  file  of  LA-ICP-MS 
U-Pb  isotopic  data  and  ages.  If  you  are  viewing  the 
PDF  of  this  paper  or  reading  it  offline,  please  visit 
http://dx.doi.org/10.1130/GES00614.Sl  or  the  full- 
text  article  on  www.gsapubs.org  to  view  Supplemen¬ 
tal  Table  1. 


zircon  data  from  other  Proterozoic  quartzite 
exposed  in  New  Mexico  in  the  statistical  analy¬ 
sis  to  evaluate  possible  correlations  throughout 
the  broader  region.  Ortega  Formation  quartzite 
exposed  in  the  Tusas  Mountains  (ORT-N;  Jones 
et  al.,  2009)  was  deposited  ca.  1690-1670  Ma 
and  derived  from  >1700  Ma  rocks  in  northern 
New  Mexico  and  southern  Colorado  (Jones  et 
al.,  2009).  Quartzite  exposed  in  the  Manzano 
Mountains  (K05-ABO-1Q;  Luther,  2006)  and 
San  Andres  Mountains  (SAHC-1Q;  Amato  et 
al.,  2008)  was  deposited  ca.  1650-1600  Ma 
and  contains  abundant  <1700  Ma  detrital  zircon 
derived  from  sources  in  northern  and  central 
New  Mexico.  Results  of  the  K-S  test  are  pre¬ 
sented  in  Table  2.  The  test  mathematically  com¬ 
pares  two  age  distributions  to  determine  if  there 
is  a  statistically  significant  difference  between 
the  two  and  returns  the  probability  (P)  that  the 
two  age  distributions  were  drawn  from  the  same 
population.  The  P  value  must  be  >0.05  to  be 
95%  confident  that  two  populations  are  not  sta¬ 
tistically  different.  The  higher  the  P  value,  the 
more  likely  it  is  that  the  two  age  distributions 
were  drawn  from  the  same  population. 

Results 

In  general,  the  detrital  zircon  morphologies 
were  remarkably  similar  among  all  samples 
and  are  dominated  by  subhedral,  subrounded 
to  rounded,  and  equant  to  slightly  elongate  or 
prismatic  grains.  A  few  notable  exceptions 
are  described  below.  Age  data  summarized  in 
Table  1  indicate  that  most  of  the  detrital  zir¬ 
con  grains  analyzed  were  Proterozoic,  with  a 
single  dominant  age  peak  between  1800  and 
1700  Ma  (Fig.  7).  Among  all  samples,  the  per¬ 
centage  of  Archean  grains  recognized  was  only 
0.9%-7.4%.  Analytical  results  described  below 
are  grouped  by  stratigraphic  association  and 
are  presented  in  order  of  previously  interpreted 
stratigraphic  succession. 

Big  Rock  Formation  and 
Marquehas  Formation 

We  collected  quartzite  and  conglomerate 
from  the  middle  and  upper  portion  of  the  Big 
Rock  and  Marquenas  Formations  to  test  the 
proposed  correlation  of  the  two  units  between 
the  Tusas  and  Picuris  Mountains.  In  the  Tusas 
Mountains,  cross-bedded  quartzite  from  the 
middle  of  the  Big  Rock  Formation  (J07-BR4q) 
was  collected  in  the  hinge  zone  of  the  Big 
Rock  syncline  (Aby  et  al.,  2010).  This  sample 
contained  the  greatest  abundance  of  euhedral, 
prismatic  grains  with  only  -5%  subrounded  to 
rounded  grains.  Detrital  zircon  “"Pb/^Pb  ages 
ranged  from  2726  to  1674  Ma,  and  the  age- 
distribution  plot  contains  a  single  broad  peak 


with  a  maximum  at  1704  Ma  together  with  a 
small  peak  at  1840  Ma  (Fig.  7H).  A  cluster  of 
three  grains  also  defines  a  small  peak  at  2666  Ma 
(Fig.  5G;  Table  1).  Cross-bedded  micaceous 
quartzite  from  the  upper  Big  Rock  Formation 
(J07-BR5)  was  collected  along  the  southeast 
limb  of  the  Posos  anticline  (Aby  et  al.,  2010). 
Detrital  zircon  2ll7Pb/206Pb  ages  from  this  sample 
range  from  3322  to  1665  Ma,  and  the  age  distri¬ 
bution  plot  has  only  a  single  broad  peak  with  a 
mode  of  1720  Ma  (Fig.  7G). 

In  the  Picuris  Mountains,  we  collected  two 
samples  of  the  Marquenas  Formation  along 
Cerro  de  las  Marquenas  between  New  Mexico 
State  Highway  75  to  the  south  and  the  Plomo 
fault  and  Copper  Hill  anticline  to  the  north  (Fig. 
3  cross  section).  Quartzite  from  the  middle 
part  of  the  unit  (J07-PIC2)  yielded  detrital  zir¬ 
con  with  207Pb/2®Pb  ages  ranging  from  2900  to 
1465  Ma,  and  the  age-distribution  curve  shows 
dominant  and  subsidiary  age  peaks  at  1707  and 
1870  Ma,  respectively  (Fig.  7B).  Quartz-pebble 
conglomerate  from  the  upper  part  of  the  unit 
(J07-PIC3)  yielded  a  detrital  zircon  207Pb/206Pb 
age  range  of  2946  to  1420  Ma  (Table  1)  with  a 
single  dominant  age  peak  at  1715  Ma  (Fig.  7A). 
However,  unlike  any  of  the  other  units  collected 
as  part  of  this  study,  the  youngest  age  peaks  in 
both  of  the  Marquenas  Formation  samples  were 
Mesoproterozoic  with  ages  of  1479  Ma  in  the 
middle  quartzite  (Fig.  7B)  and  1457  Ma  in  the 
upper  conglomerate  (Fig.  7A).  The  Mesopro¬ 
terozoic  peaks  are  defined  by  six  out  of  222 
grains  in  the  middle  quartzite  (2.7%,  Table  1) 
and  21  out  of  199  grains  in  the  upper  conglom¬ 
erate  (10.6%,  Table  1). 

The  Paleoproterozoic  portions  of  the  age- 
distribution  spectra  for  the  Big  Rock  and  Mar¬ 
quenas  Formations  overlap  with  relatively  minor 
variation  among  the  peak  ages.  The  main  differ¬ 
ences  between  the  four  curves  are  the  broadness 
of  the  two  Big  Rock  Formation  peaks  compared 
to  the  relatively  narrow  Marquenas  Formation. 
The  dominant  peaks  in  the  Marquenas  Formation 
samples  also  skew  toward  slightly  younger  ages, 
indicating  a  greater  abundance  of  <1700  Ma 
grains.  Results  of  the  statistical  tests  shown  in 
Table  2  indicate  that  the  two  units  of  the  Mar¬ 
quenas  Formation  are  significantly  similar  with 
a  P  value  of  0.357.  However,  only  the  Upper 
Marquenas  Formation  is  similar  to  both  units 
of  the  Big  Rock  Formation  (P  values  of  0.100 
and  0.179).  The  Middle  Marquenas  Formation 
is  statistically  similar  to  the  Middle  Big  Rock 
Formation  (P  =  0.242)  but  not  the  upper  part  of 
the  unit  (P  =  0.034). 

The  discovery  of  ca.  1479-1459  Ma  zircon 
in  the  two  Marquenas  Formation  samples  was 
the  most  surprising  result  of  this  study.  We  col¬ 
lected  backscattered  electron  (BSE)  images 
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Figure  5.  Age-distribution  diagrams  for  2ll7Pb/20fiPb  ages  and  U-Pb  concordia  diagrams  for  detrital  zircon  from  the  Marque- 
nas  and  Big  Rock  formations,  Picuris  and  Tusas  Mountains,  New  Mexico.  Plots  represent  all  grains  analyzed,  and  ages  of 
significant  peaks  are  indicated  (Ma,  Table  1).  Histogram  bins  are  25  m.y.  wide.  See  Table  1  for  a  summary  of  age  data  and 
Supplemental  Table  1  (see  footnote  1)  for  complete  data  set. 
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of  the  Mesoproterozoic  grains  to  determine 
if  the  young  ages  were  the  result  of  metamor- 
phic  overgrowths,  but  most  of  the  grains  had 
well-developed  concentric  growth  zones  and  no 
noticeable  overgrowths  or  other  internal  irregu¬ 
larities  (Fig.  8).  The  other  grains  did  not  produce 
any  recognizable  internal  patterns.  The  subhe- 
dral  external  morphology  of  the  ca.  1450  Ma 
grains  was  consistent  with  other  detrital  zircon 
grains  from  the  same  sample  and  from  nearly  all 
of  the  other  samples  analyzed  during  this  study. 
The  grains  do  not  match  the  external  morphology 
of  typical  metamorphic  zircon  shown  in  Corfu  et 
al.  (2003).  The  Th/U  ratios  of  the  ca.  1450  Ma 
grains  were  among  the  lowest  measured  for  the 
two  samples  with  a  range  of  1.24—0.24  ( Supple¬ 
mental  Table  1  [see  footnote  1]),  but  they  are  not 


significantly  different  from  the  rest  of  the  pop¬ 
ulation  or  other  samples  and  do  not  match  the 
very  low  ratios  that  are  typical  for  metamorphic 
zircon  (<0.07;  Rubatto,  2002).  When  data  acqui¬ 
sition  for  the  two  Marquenas  Formation  samples 
was  initially  completed  in  2009,  we  reanalyzed 
a  subset  of  the  ca.  1450  Ma  grains  in  the  upper 
conglomerate  to  confirm  the  original  age  deter¬ 
minations.  In  all  but  one  case,  both  ages  were 
identical  within  analytical  uncertainty  (Figs. 
8A-8C).  In  the  other  case,  the  second  analysis 
produced  an  Archean  age,  suggesting  that  the 
grain  had  an  inherited  component  in  its  core. 
Only  one  ca.  1470  Ma  detrital  zircon  grain  was 
initially  recognized  in  the  middle  Marquenas 
Formation  sample.  We  picked  and  remounted 
a  new  set  of  zircon  samples  from  the  two  Mar¬ 


quenas  Formation  samples  that  were  analyzed  in 
2010  using  the  same  techniques,  and,  again,  both 
samples  yielded  ages  <1500  Ma.  The  proportion 
of  ca.  1470-1450  Ma  grains  was  -10%  of  the 
total  in  both  data  sets  for  the  upper  conglomer¬ 
ate  (  Supplemental  Table  1  [see  footnote  1]),  and 
the  second  set  of  analyses  yielded  a  statistically 
significant  population  of  ca.  1470  Ma  grains  in 
the  middle  quartzite  (four  out  of  115  analyzed; 
Supplemental  Table  1  [see  footnote  1]). 

Ortega  Formation 

We  collected  two  samples  from  the  Ortega 
Formation  in  the  Picuris  Mountains  to  test  the 
correlation  with  exposures  in  the  Tusas  Moun¬ 
tains  and  to  evaluate  possible  changes  in  detrital 
zircon  populations  along  strike  and  vertically 
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Figure  6.  Age-distribution  diagrams  for  “’Pb/^Pb  ages  and  U-Pb  concordia  diagrams  for  detrital  zircon  from  the 
Ortega  and  Rinconada  formations  of  the  Hondo  Group,  Picuris  Mountains,  New  Mexico.  Plots  represent  all  grains 
analyzed,  and  ages  of  significant  peaks  are  indicated  (Ma,  Table  1).  Histogram  bins  are  25  m.y.  wide.  See  Table  1  for 
a  summary  of  age  data  and  Supplemental  Table  1  (see  footnote  1)  for  complete  analytical  data. 
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within  the  stratigraphic  succession.  One  sample 
of  Ortega  Formation  quartzite  (J07-PIC1),  col¬ 
lected  from  exposures  in  the  northern  limb  of 
the  Hondo  Syncline  (Fig.  3),  yielded  detrital  zir¬ 
con  with  “PbPPb  ages  ranging  from  2686  to 
1670  Ma.  The  dominant  age  peak  was  1724  Ma 
(Fig.  7F).  Three  grains  define  a  separate  small 
age  peak  at  2483  Ma  (Fig.  6E;  Table  1).  Quartz¬ 
ite  from  the  uppermost  part  of  the  unit  (J07- 
PIC4)  was  collected  in  the  core  of  the  Cop¬ 
per  Hill  anticline  (Fig.  3  cross  section),  and  it 
yielded  detrital  zircon  with  207Pb/206Pb  ages  of 
2658-1538  Ma  (Table  1).  The  age  spectrum 
shows  a  single  broad  peak  with  maxima  at  1729 
and  1765  Ma  (Fig.  7D). 

The  age-distribution  curve  and  peak  age  for 
detrital  zircon  from  the  middle  Ortega  Forma¬ 
tion  in  both  the  Picuris  and  Tusas  Mountains  are 
nearly  identical  (Figs.  7E  and  7F).  The  statisti¬ 
cal  test  also  indicates  that  these  two  samples 
are  significantly  similar  with  a  relatively  high 
P  value  of  0.425  (Table  2).  The  middle  Ortega 
Formation  in  the  Picuris  Mountains  is  also 
statistically  similar  with  the  middle  Big  Rock 
Formation  (P  =  0.523),  but  P<0.05  when  com¬ 
pared  with  all  other  samples  (Table  2).  The 
upper  Ortega  Formation  is  significantly  simi¬ 
lar  with  the  upper  Marquenas  Formation  (P  = 
0.351)  and  all  three  samples  of  Ortega  and  Big 
Rock  Formation  from  the  Tusas  Mountains  (P  = 
0.999-0.596),  with  a  P  value  of  nearly  1  when 


Figure  7.  Summary  of  detrital  zircon 
207Pb/2#6Pb  ages  (<10%  discordant,  Protero¬ 
zoic  grains)  and  basement  U-Pb  crystalliza¬ 
tion  ages  in  Proterozoic  exposures  in  north- 
central  New  Mexico.  Relative  probability 
diagrams  for:  (A)  conglomerate  of  the  upper 
Marquenas  Formation,  Picuris  Mountains; 
(B)  quartzite  of  the  middle  Marquenas  For¬ 
mation,  Picuris  Mountains;  (C)  quartzite 
of  the  upper  Rinconada  Formation,  Picuris 
Mountains;  (D)  quartzite  of  the  upper 
Ortega  Formation,  Picuris  Mountains; 
(E)  quartzite  of  the  middle  Ortega  Forma¬ 
tion,  Tusas  Mountains,  previously  published 
by  Jones  et  al.  (2009);  (F)  quartzite  of  the 
middle  Ortega  Formation,  Picuris  Moun¬ 
tains;  (G)  quartzite  of  the  upper  Big  Rock 
Formation,  Tusas  Mountains;  (H)  quartzite 
of  the  middle  Big  Fork  Formation,  Tusas 
Mountains;  (I)  basement  U-Pb  crystalliza¬ 
tion  ages  for  exposures  in  northern  and 
central  New  Mexico  (Karlstrom  et  al.,  2004, 
and  references  therein).  Ages  of  significant 
peaks  are  indicated  (Ma,  Table  1).  Curves 
are  normalized  so  that  area  underneath 
each  is  the  same  (Gehrels,  2009). 
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J07-PIC2 

J07-PIC5 

J07-PIC4 

ORT-N 

J07-PIC1 

J07-BR5 

J07-BR4q 

K05- ABO-1  Q 

SAHC-1Q 


TABLE  2.  DETRITAL  ZIRCON  AGE  SPECTRA  COMPARISON  USING  OVERLAP,  SIMILARITY,  AND  K-S  STATISTICS 

_ Overlap _ 

J07-PIC3  J07-PIC2  J07-PIC5  J07-PIC4  ORT-N  J07-PIC1  J07-BR5  J07-BR4Q 


0.722 

0.728  0.787 

0.847  0.593 

0.448  0.644 

0.595  0.642 

0.626  0.600 

0.545  0.693 

0.599  0.738 

0.362  0.483 


0.610 

0.586  0.452 

0.794  0.555 

0.648  0.584 

0.709  0.465 

CL509  CL424 


0.677 

0.557  0.653 

0.774  0.707 

0.646  0.706 

0.588  0.551 


0.546 

0.577  0.696 

0.416  0.530 


K05-ABO-1Q 


0.643 


J07-PIC2 

J07-PIC5 

J07-PIC4 

ORT-N 

J07-PIC1 

J07-BR5 

J07-BR4q 

K05- ABO-1  Q 

SAHC-1Q 


_ Similarity _ 

J07-PIC3  J07-PIC2  J07-PIC5  J07-PIC4  ORT-N  J07-PIC1  J07-BR5 

0.897 

0.833  0.842 


0.838 

0.820 

0.862 

0.842 

0.819 

0.774 


0.882 

0.855 

0.891 

0.897 

0.908 

0.872 


0.870 

0.793 

0.868 

0.863 

0.756 

0.678 


0.907 

0.840  0.893 

0.916  0.903  0.857 

0.888  0.915  0.885 

0.858  0.805  0.804 

0.824  0.754  0.753 


0.887 


(X801 


0.841 

0.781 


K05-ABO-1Q 


0.928 


J07-PIC3  J07-PIC2  J07-PIC5 

jo7-pic2 

J07-PIC5  0.000  0.000 

J07-PIC4  0.019  0.039 

ORT-N  0.038  0.009  0.000 

J07-PIC1  0.030  0.014  0.000 

J07-BR5  mill  0.017 

J07-BR4Q  0.100  0.242  0.000 

K05- ABO-1  Q  0.004  0.000 

SAHC-1Q_ 0001_ 0.000_ 0.000 


K-S  statistic _ 

J07-PIC4  ORT-N  J07-PIC1  J07-BR5 


0.012  0.425 

0.999  0.890  0.031 

0.596  0.756  0.523  0.835 

0.000  0.000  0.000  0.000 

0.000  0.000_ 0.000_ 0.000 


0.001  _ 

0.000_ 0.624 


Note :  Overlap  numbers  indicate  the  degree  of  overlap  of  two  age  probabilities  where  1 .0  is  perfect  overlap  and  0.0  is  no  overlap.  Similarity  numbers 
are  a  measure  of  similarity  in  the  proportion  of  overlapping  ages  where  higher  values  up  to  1 .0  reflect  similar  proportions  of  overlapping  ages  and  lower 
values  down  to  0.0  reflect  different  proportions  of  ages  that  may  or  may  not  overlap.  K-S  statistic  numbers  indicate  the  probability  (P)  that  two  samples 
are  derived  from  the  same  population.  The  higher  the  value,  the  more  likely  it  is  that  the  two  age  distributions  were  drawn  from  the  same  population.  The 
P  value  must  exceed  0.05  to  be  95%  confident  that  the  two  populations  are  not  statistically  different,  and  those  values  exceeding  0.05  are  shown  in  bold 
and  shaded  in  gray.  Statistics  were  calculated  using  Microsoft  Excel  macros  made  available  by  G.  Gehrels  at  the  University  of  Arizona  Laserchron  Center 
(Gehrels,  2009).  Samples  ORT-N  (Jones  et  al.,  2009),  K05-ABO-1Q  (Luther,  2006),  and  SAHC-1Q  (Amato  et  al.,  2008)  are  Paleoproterozoic  quartzites 
exposed  in  the  Tusas  Mountains,  Manzano  Mountains,  and  San  Andres  Mountains,  respectively,  and  are  included  for  comparison  and  discussed  in  the 
accompanying  text. 


compared  with  the  upper  Big  Rock  Formation. 
However,  P<0.05  when  compared  with  all  other 
samples  (Table  2). 

Rinconada  Formation 

We  collected  one  sample  (J07-PIC5)  from 
the  upper  quartzite  (Bauer  and  Helper,  1994) 
of  the  Rinconada  Formation  (Hr5)  exposed  on 
the  south  limb  of  the  Copper  Hill  anticline  (Fig. 
3  cross  section)  to  evaluate  possible  changes  in 
the  detrital  zircon  population  vertically  within 
the  Hondo  Group  stratigraphic  succession  and 
to  provide  a  baseline  for  correlation  of  Rincon¬ 
ada  Formation  units  in  exposures  throughout  the 
surrounding  region  (Fig.  2).  This  sample  con¬ 
tained  the  greatest  abundance  of  rounded  to  sub¬ 
rounded,  equant  zircon  and  few,  if  any,  elongate 
grains.  Detrital  zircon  207Pb/206Pb  ages  ranged 
from  2892  to  1686  Ma  (Table  1),  and  the  age- 
distribution  curve  is  dominated  by  a  single  peak 
with  a  maximum  age  of  1762  Ma  (Fig.  7C).  A 
cluster  of  five  Archean  grains  also  defines  a 
small  peak  at  2592  Ma  (Fig.  6A;  Table  1). 

The  age-distribution  curve  for  the  Upper 
Rinconada  Formation  is  notably  different  from 


the  other  samples  in  that  it  lacks  a  peak  at  1729- 
1704  Ma  and  skews  toward  slightly  older  ages. 
The  age  distribution  is  most  similar  to  the  Upper 
Ortega  Formation  and  both  Big  Rock  Formation 
samples  because  all  three  have  a  ca.  1760  Ma 
component.  However,  the  Upper  Rinconada 
Formation  is  not  statistically  similar  with  any  of 
the  other  samples  included  in  this  study.  Com¬ 
parison  with  two  samples,  upper  Ortega  Forma¬ 
tion  and  upper  Big  Rock  Formation,  produces  P 
values  of  0.039  and  0.017,  respectively,  whereas 
comparison  with  all  other  samples  returns  val¬ 
ues  of  zero  (Table  2). 

DISCUSSION 

Provenance  Trends  of  the  Yadito  and 
Hondo  Groups 

Our  results  indicate  that  the  detrital  zircon 
populations  are  very  similar  among  all  samples 
of  the  Vadito  and  Hondo  groups.  We  interpret 
the  relatively  small  number  of  Archean  grains  to 
represent  detrital  zircon  derived  from  the  Wyo¬ 
ming  Province  to  the  north  (Chamberlain  et  al.. 


2003;  Jones  et  al.,  2009),  but  they  might  rep¬ 
resent  reworked  Paleoproterozoic  sedimentary 
rocks  (Hill  and  Bickford,  2001;  Premo  et  al., 
2007).  The  1876-1840  Ma  peaks  that  are  appar¬ 
ent  in  only  three  samples  might  represent  zircon 
derived  from  the  Trans-Hudson  orogen  or  Black 
Hills  to  the  north  (Van  Schmus  et  al.,  1987; 
Redden  et  al.,  1990)  or  from  the  Grand  Canyon 
region  or  Mojave  region  to  the  west  (Wooden  et 
al.,  1988;  Wooden  and  DeWitt,  1991;  Hawkins 
et  al.,  1996;  Shufeldt  et  al.,  2010).  Otherwise, 
the  dominant  peaks  on  the  age-distribution 
curves  overlap  with  zircon  crystallization  ages 
of  basement  exposures  from  throughout  the  sur¬ 
rounding  region  (Fig.  71).  Older  peak  ages  of 
ca.  1765-1730  Ma  represent  zircons  that  were 
derived  from  plutonic  complexes  or  greenstone 
associations  such  as  the  Pecos,  Moppin,  or  Gold 
Hill  successions  in  northern  New  Mexico  (Fig. 
2;  Karlstrom  et  al.,  2004,  and  references  therein). 
Younger  peak  ages  of  ca.  1730-1704  Ma 
were  likely  derived  from  volcano-plutonic  suc¬ 
cessions  that  are  exposed  throughout  south¬ 
ern  Colorado  and  northern  New  Mexico  (e.g., 
Bickford  et  al.,  1989).  We  interpret  these 
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Figure  8.  Backscattered  electron  (BSE)  images  (A-D)  for  selected  ca.  1470-1450  Ma 
grains  from  the  upper  Marqueiias  Formation  conglomerate  showing  laser  ablation 
pits  and  associated  ages  (Ma).  The  photo  mosaic  on  the  right  (E)  shows  a  representa¬ 
tive  portion  of  the  grain  mount  for  sample  J07-PIC3  that  contains  multiple  <1500  Ma 
grains.  All  grain  numbers  represent  the  analysis  number  and  are  indexed  to  the  data 
presented  in  Supplemental  Table  1  (see  footnote  1). 


relationships  to  indicate  predominately  local 
sources  for  all  of  the  units,  with  variations  in  age 
distributions  and  age  peaks  corresponding  to  the 
age  of  the  underlying  basement.  Our  findings 
are  consistent  with  Pb  isotopic  data  of  McLen¬ 
nan  et  al.  (1995)  and  limited  detrital  zircon  data 
from  the  Ortega  Formation  reported  by  Aleini- 
koff  et  al.  (1993).  Our  findings  also  expand  the 
detrital  zircon  data  set  presented  by  Jones  et  al. 
(2009)  for  ca.  1.70  Ga  metasedimentary  succes¬ 
sions  formed  in  the  southern  Yavapai  Province 
and  support  that  model  for  relatively  short-lived 
basins  that  were  dominated  by  locally  derived 
Paleoproterozoic  sediment. 

Similarities  in  the  detrital  zircon  age  spec¬ 
tra  for  samples  from  the  Tusas  Mountains  are 
consistent  with  local  evidence  for  a  gradational 
depositional  contact  between  the  Vadito  and 
Hondo  groups  (Bauer  and  Williams,  1989).  The 
age  ranges  and  peak  ages  are  similar  among 
the  three  samples,  though  the  peak  in  the  age- 
distribution  curve  is  narrower  upsection  in  the 
Ortega  Formation  quartzite  than  in  the  Big  Rock 
Formation  below  (Fig.  7).  This  trend  might  be 
the  result  of  changing  local  source  rocks,  with 
highly  variable  initial  sources  giving  way  to 
more  homogeneous  sources  as  the  depositional 
systems  evolved.  The  change  might  also  reflect 
contrasting  depositional  environments,  with  the 
thick  pebble  to  boulder  conglomerate  of  the  Big 
Rock  Formation  indicating  high-energy  fluvial 
systems  and  the  massive  quartzite  of  the  Ortega 
Formation  representing  a  shallow  marine  envi¬ 
ronment  (Soegaard  and  Eriksson,  1985).  The 
age-distribution  curves  for  middle  Ortega  For¬ 
mation  quartzite  are  nearly  identical  between  the 
Tusas  and  Picuris  Mountains,  supporting  the  cor¬ 
relation  of  Hondo  Group  rocks  between  the  two 
ranges.  The  shift  in  peak  ages  toward  ca.  1760  Ma 
in  the  upper  Ortega  and  Rinconada  Formations 
indicates  another  change  in  source  rocks  toward 
the  end  of  Ortega  Formation  deposition. 

Age  and  Provenance  of 
the  Marqueiias  Formation 

Depositional  Age  Constraints 

We  interpret  the  results  described  above  to 
indicate  that  the  two  samples  from  the  Mar- 
quenas  Formation  have  ca.  1479-1457  Ma 
detrital  zircon,  a  finding  that  has  never  been 
reported  among  presumed  Paleoproterozoic 
metasedimentary  rocks  throughout  the  region. 
The  2(l7Pb/2,,6Pb  ages  of  the  five  youngest  grains 
from  the  middle  quartzite  range  from  1489  to 
1465  Ma,  with  a  weighted  mean  of  1477  ± 
13  Ma.  The  207Pb/206Pb  ages  of  23  grains  from 
the  upper  conglomerate  range  from  1481  to 
1420  Ma,  with  a  weighted  mean  of  1453  ± 
10  Ma.  These  ages  are  interpreted  as  the  best 
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estimates  for  the  maximum  depositional  age  of 
the  middle  and  upper  Marquenas  Formation, 
respectively.  These  new  ages  indicate  that  the 
Marquenas  Formation  is  not  correlative  with 
the  ca.  1700  Ma  Vadito  Group,  and  our  reinter¬ 
preted  stratigraphic  relationships  are  shown  in 
Figure  4C.  The  northern  contact  of  the  Marque¬ 
nas  Formation  is  the  Plomo  fault,  a  well-defined 
ductile  shear  zone.  The  southern  margin  is  not 
well  exposed  but  must  represent  an  unconfor¬ 
mity,  fault,  or  both.  It  is  possible  that  other  rocks 
included  within  the  Vadito  Group  are  associ¬ 
ated  with  the  younger  Marquenas  Formation 
as  well;  however,  more  work  is  needed  to  test 
this  possibility.  The  minimum  depositional  age 
of  the  Marquenas  Formation  is  not  constrained 
directly  by  cross-cutting  relationships.  The 
Picuris  Pueblo  granite  exposed  in  the  southeast¬ 
ern  Picuris  Mountains  (Bauer  et  al.,  2005)  does 
appear  to  intrude  eastern  exposures  of  the  Mar¬ 
quenas  Formation  (Fig.  3)  and  may  provide  an 
independent  age  constraint  on  deposition,  but  it 
has  not  yet  been  dated. 

The  best  available  constraints  on  the  depo¬ 
sitional  age  of  the  Marquenas  Formation  come 
from  previous  detailed  structural  and  metamor- 
phic  studies  of  the  Picuris  and  adjacent  moun¬ 
tain  ranges  that  are  described  above.  All  stud¬ 
ies  of  the  Picuris  Mountains  conclude  that  the 
Marquenas  Formation  experienced  the  same 
regional-scale  F2  folding  event  as  adjacent 
rocks  of  the  Vadito  and  Hondo  groups  (Miller 
et  al.  1963;  Nielsen  and  Scott,  1979;  Holcombe 
and  Callender,  1982;  Mawer  et  al.  1990;  Bauer, 
1993;  Williams  et  al.  1999).  Furthermore, 
there  is  abundant  evidence  that  F2  folding  was 
broadly  contemporaneous  with  amphibolite- 
facies  metamorphism  involving  the  growth  of 
all  three  aluminosilicate  minerals,  and  metamor- 
phic  minerals  are  locally  present  in  aluminum- 
rich  layers  of  the  Marquenas  Formation.  Thus, 
the  metamorphic  ages  described  above  provide 
a  reasonable  lower  age  constraint  on  the  Mar¬ 
quenas  Formation  and  suggest  that  the  unit  was 
deposited,  deformed,  and  metamorphosed  dur¬ 
ing  the  interval  between  ca.  1453  and  1435  Ma. 

Provenance  of  the  Marquenas  Formation 

A  ca.  1453  Ma  or  younger  depositional  age 
means  that  the  Marquenas  Formation  could 
be  derived  from  the  Vadito  and  Hondo  groups, 
consistent  with  the  interpretations  of  Soegaard 
and  Eriksson  (1986).  However,  Bauer  and  Wil¬ 
liams  (1989)  and  Mawer  et  al.  (1990)  pointed 
out  the  disparity  of  aluminum  silicate  miner¬ 
als  in  Marquenas  Formation  cobbles  relative  to 
Ortega  Formation  quartzite  and  argued  that  the 
Ortega  Formation  was  not  a  likely  source  for 
the  quartzite  clasts.  Derivation  of  detrital  zir¬ 
con  from  the  Vadito  and  Hondo  groups  would 


help  to  explain  the  strong  similarity  of  detrital 
zircon  populations  between  the  upper  Marque¬ 
sas  Formation  and  all  other  samples  except 
for  the  upper  Rinconada  Formation.  However, 
both  Marquenas  Formation  samples  have  abun¬ 
dant  <1700  Ma  zircon  (Fig.  7;  Table  1)  that 
are  not  present  in  the  Vadito  or  Hondo  Group. 
We  interpret  these  ages  to  represent  ca.  1700- 
1650  Ma  volcano-plutonic  successions  that 
dominate  Proterozoic  exposures  in  northern 
and  central  New  Mexico  and  include  multiple 
granitic  plutons  that  intrude  the  Vadito  Group  in 
the  southern  Picuris  Mountains  (Figs.  2  and  3; 
Karlstrom  et  al.,  2004,  and  references  therein). 

Similar  detrital  zircon  ages  are  reported  for 
ca.  1650-1600  Ma  quartzites  exposed  to  the 
south  in  the  Manzano  and  San  Andres  Moun¬ 
tains  (Luther,  2006;  Amato  et  al.,  2008),  and  the 
K-S  test  confirms  the  similarity  of  detrital  zircon 
populations  among  these  exposures  (P  =  0.624; 
Table  2)  and  between  the  middle  Marquenas 
Formation  and  Manzano  Mountains  exposures 
(P  =  0.054;  Table  2).  The  dissimilarity  between 
the  upper  Marquenas  Formation  and  Manzano 
Mountains  quartzite  (P  =  0.004)  could  indicate 
a  shift  to  more  local  or  northern  >1700  Ma 
sources  upsection,  an  interpretation  that  is  sup¬ 
ported  by  the  slightly  older  peak  age  in  the  upper 
conglomerate  and  more  pronounced  asymmetry 
on  the  right-hand  (i.e.,  older)  side  of  the  domi¬ 
nant  age  peak  (Fig.  7A). 

The  range  of  Mesoproterozoic  ages  in  the 
two  Marquenas  Formation  samples  corresponds 
well  with  the  ages  of  granitoids  that  intruded 
throughout  the  region  1470-1360  Ma  (Fig.  71; 
Reed  et  al.,  1993;  Karlstrom  et  al.,  2004)  and 
overlaps  with  widespread  1480-1380  Ma  rhy¬ 
olites  in  the  subsurface  of  the  midcontinent  to 
the  east  (Van  Schmus  and  Bickford,  1981).  The 
1479  and  1457  Ma  age  peaks  for  the  middle 
and  upper  Marquenas  Formation  overlap  the 
upper  end  of  the  age  spectrum  for  Mesoprotero¬ 
zoic  igneous  rocks  but  predate  the  ca.  1440- 
1430  Ma  magmatic  peak  in  the  Rocky  Moun¬ 
tains  and  southwestern  United  States  (Reed  et 
al.,  1993;  Karlstrom  et  al.,  2004). 

Regional  Tectonic  Implications 

The  discovery  of  Mesoproterozoic  detrital 
zircon  in  the  presumed  Paleoproterozoic  Mar¬ 
quenas  Formation  is  significant  in  a  regional 
tectonic  context  because  it  requires  that  D2 
deformation  in  northern  New  Mexico  occurred 
ca.  1 .43  Ga  rather  than  ca.  1 .65  Ga  as  previously 
interpreted.  We  contend  that  the  Marquenas  For¬ 
mation  represents  the  first  known  occurrence 
of  pre-  to  synorogenic  metasedimentary  rocks 
associated  with  ca.  1.4  Ga  intracontinental  oro¬ 
genesis  in  the  southwestern  United  States.  The 


Marquenas  Formation  is  truncated  to  the  west 
along  the  Plomo  fault  but  may  be  equivalent  to 
conglomerate  exposed  along  strike  to  the  east 
for  -60  km  in  the  Rio  Mora  region  (Grambling 
and  Codding,  1982).  These  patterns  suggest  that 
deposition  could  have  been  locally  controlled 
by  movement  along  the  Plomo  fault  or  a  pre¬ 
cursor  structure.  The  Marquenas  Formation 
was  likely  deposited  unconformably  on  Vadito 
Group  exposures,  thus  explaining  the  apparent 
stratigraphic  succession  south  of  the  Plomo  fault 
(Bauer,  1993).  Alternatively,  the  Marquenas 
Formation  was  faulted  against  the  Vadito  Group 
prior  to  both  rock  units  being  folded  and  faulted 
against  the  Hondo  Group.  Either  interpretation 
requires  significant  deformation  of  the  Marque¬ 
sas  Formation  after  deposition  and  provides 
additional  evidence  that  much  of  the  deforma¬ 
tion  and  metamorphism  observed  in  the  Picuris 
Mountains  occurred  ca.  1.4  Ga.  Furthermore, 
our  findings  provide  a  new  line  of  evidence  in 
the  case  for  ca.  1.4  Ga  regional  orogenesis  in 
the  southern  Rocky  Mountains  within  a  conver¬ 
gent  or  transpressional  tectonic  setting  (Nyman 
et  al.,  1994;  Kirby  et  al.,  1995;  Selverstone  et 
al.,  2000;  Shaw  et  al.,  2001;  Jessup  et  al.,  2005; 
Jones  et  al.,  2010). 

We  envision  a  scenario  in  which  the  boulder 
conglomerate  and  associated  quartzite  were 
deposited  in  a  broad  alluvial  fan  that  formed  in 
response  to  a  growing  orogenic  highland  to  the 
south.  These  highlands  may  be  associated  with 
ca.  1.4  Ga  movement  along  the  north- vergent 
Manzano  thrust  belt  (Rogers,  2001;  Baer  et 
al.,  2003),  a  series  of  northeast-striking,  north- 
vergent  thrust  faults  and  ductile  shear  zones 
exposed  in  central  New  Mexico  (Fig.  1;  Karl¬ 
strom  et  al.,  2004;  Cather  et  al.,  2006).  Uplift 
of  the  Manzano  Group,  a  thick  succession  of 
volcanic  and  sedimentary  rocks  (Luther  et  al., 
2005b;  Luther,  2006),  along  the  thrust  belt  pro¬ 
vided  ample  sources  of  Paleoproterozoic  detrital 
zircon  as  well  as  abundant  quartzite  and  rhyolite 
clasts.  Alternatively,  the  Marquenas  Formation 
was  derived  from  more  local  sources  such  as 
the  Pecos  greenstone,  lower  Vadito  Group,  or 
the  many  ca.  1670-1620  Ma  and  ca.  1450  Ma 
granites  exposed  in  the  surrounding  region  (Fig. 
2;  Karlstrom  et  al.,  2004).  In  this  scenario,  the 
Marquenas  Formation  would  be  analogous  to 
younger  Mesoproterozoic  successions  such  as 
the  Hazel  Formation,  an  orogenic  clastic  wedge 
deposited  ca.  1 120-1080  Ma  during  the  Gren¬ 
ville  orogeny  in  the  Van  Horn  area  of  western 
Texas  (Fig.  1;  Soegaard  and  Callahan,  1994; 
Mosher,  1998;  Bickford  et  al.,  2000).  The  Hazel 
Formation  is  made  up  of  2500  m  of  immature 
boulder  conglomerate  and  fine-grained  sand¬ 
stone  that  were  derived  from  a  southern  tectonic 
highland  and  deposited  in  large,  coalescing 
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alluvial  fans  immediately  north  of  the  frontal 
zone  of  the  orogen  (Soegaard  and  Callahan, 
1994).  Age  constraints  indicate  that  the  boulder 
conglomerate  was  deposited  within  a  ca.  50  m.y. 
time  window  prior  to  deformation,  a  time  frame 
that  is  similar  to  the  constraints  for  the  Marque¬ 
sas  Formation  described  above. 

CONCLUSIONS 

Detrital  zircon  geochronology  provides  a  crit¬ 
ical  test  of  correlations  of  Proterozoic  metasedi¬ 
mentary  rocks  exposed  within  the  >1000  km2 
Pilar  basin  in  northern  New  Mexico  and  yields 
new  insight  into  their  sedimentary  provenance. 
Our  results  support  previous  correlations  of 
quartzite  of  the  Hondo  Group  between  the 
Tusas  and  Picuris  Mountains  and  indicate  that 
they  were  predominately  derived  from  local 
Paleoproterozoic  sources.  Systematic  shifts  in 
the  peak  ages  of  the  detrital  zircon  population 
within  the  Hondo  Group  are  interpreted  to  rep¬ 
resent  changing  source  terranes,  but  the  overall 
age  range  of  the  dominant  detrital  zircon  popula¬ 
tion  closely  matches  the  crystallization  ages  of 
basement  units  in  the  surrounding  region.  Our 
results  also  require  a  revision  of  the  Proterozoic 
lithostratigraphy  of  northern  New  Mexico.  Boul¬ 
der  conglomerate  and  quartzite  of  the  Marque¬ 
sas  Formation  were  previously  correlated  with 
the  1700  Ma  Vadito  Group,  but  detrital  zircon 
ages  require  that  it  was  deposited  after  1453  Ma. 
We  interpret  the  Marquenas  Formation  to  repre¬ 
sent  a  synorogenic  succession  that  was  derived 
from  the  south  and  deposited  ca.  1450-1435  Ma 
contemporaneous  with  regional  deformation, 
metamorphism,  and  magmatism  in  the  southern 
Rocky  Mountains.  Deformation  and  metamor¬ 
phism  of  the  Marquenas  Formation  involved  jux¬ 
taposition  against  poly-deformed,  amphibolite- 
facies  rocks  of  the  Hondo  Group  along  the  duc¬ 
tile  Plomo  fault  and  must  have  occurred  shortly 
after  deposition.  This  finding  represents  a  new 
line  of  evidence  in  support  of  regional  arguments 
for  orogenesis  at  ca.  1.4  Ga  in  a  convergent  or 
transpressional  tectonic  setting. 
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